Observation of temperature peaks due to strong viscous heating in a dusty plasma flow 
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Profound temperature peaks are observed in regions of high velocity shear in a 2D dusty plasma 
experiment with laser-driven flow. These are attributed to viscous heating, which occurs due to colli- 
sional scattering in a shear flow. Using measurements of viscosity, thermal conductivity, and spatial 
profiles of flow velocity and temperature, we determine three dimensionless numbers: Brinkman Br 
= 0.5, Prandtl Pr = 0.09, and Eckert Ec = 5.7. The large value of Br indicates significant viscous 
heating that is consistent with the observed temperature peaks. 

PACS numbers: 52.27.Lw, 52.27.Gr, 44.10. +i, 83.50.Ax 



Viscous heating occurs when there is a shear in the 
velocity profile in a fluid. Viscous heating is the energy 
lost whenever there is a gradient in flow velocity, for ex- 
ample due to moving solid surfaces such as pistons in an 
oil-lubricated engine, a human swimming in water, or an 
airplane flying through air. Despite the importance of 
this energy loss, viscous heating often does not result in 
a measurable increase in temperature, because most flu- 
ids conduct heat rapidly away from the location of flow 
gradients where the heat is generated. This competition 
between viscous heating and thermal conduction is mea- 
sured by the dimensionless parameter of fluid dynamics 
called the Brinkman number [TJ [5] , 
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Here, rj is the shear viscosity and k is the thermal conduc- 
tivity [3 , while Aw and AT are differences in flow velocity 
and temperature between two positions, for example be- 
tween a boundary and the hottest point within the flow. 
Note that Br depends on two instrinsic properties of the 
substance, viscosity 77 and thermal conductivity /v, (IS well 
as properties of the flow profile Av and AT. A substan- 
tial temperature increase due to viscous heating would 
be indicated by a value of Br of order unity [2] . While a 
value as large as unity is theoretically achievable in some 
simple flows, such a large value is seldom measured ex- 
perimentally. Most commonly, Br is orders of magnitude 
smaller than unity [H [5] . 

The Brinkman number can also be expressed as the 
product Br = PrEc of two other dimensionless numbers: 
the Prandtl number Pr and Eckert number [5] Ec, 
which are defined as 
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Here, v = rj/p is the momentum diffusivity (i.e., kine- 
matic viscosity), x = K /i c P) is the thermal diffusivity, 



c is the specific heat, and p is the mass density. The 
Prandtl number is essentially a ratio of momentum diffu- 
sivity to energy diffusivity 0], while the Eckert number 
characterizes dissipation as the ratio of a directed kinetic 
energy to a thermal kinetic energy [6] . 

The meaning of the Brinkman number is motivated by 
the energy continuity equation [7J [5] for steady condi- 
tions, 
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written here for the simplified geometry of a Couette flow. 
A Couette flow [9j is a laminar flow between two planar 
boundaries moving in the x direction at different speeds 
Av. The Brinkman number reflects the ratio of viscous 
heating to thermal conduction, i.e., the ratio of the first 
two terms in Eq. Q. The last term P cxt represents any 
external source or sink of heat energy per unit mass. 

Although the Brinkman number is usually tiny, we will 
show that it has a large value for the substance used in 
our experiment: a dusty plasma. A dusty plasma is a 
four-component mixture of micron-size particles of solid 
matter, electrons, positive ions, and a rarefied gas of neu- 
tral atoms [TOl - fM] . The small particles, which we term 
dust particles, gain a high negative electric charge Q, 
so that they collide with one another electrically, like 
charged colloids in an aqueous suspension [15] . The 
charged dust particles can self-organize into a lattice with 
crystalline properties, which can be melted to make a 
liquid [ToTl22j . A dusty plasma can be manipulated by 
external forces, which can account for P cxt in Eq. Q, to 
cause a laminar flow of dust particles [2"3"j . 

Compared with other substances, dusty plasma has 
several extreme properties. The typical interparticle 
spacing is much larger than the particle size, about 
0.5 mm as compared to 8.09 pm for our experiment, and 
correspondingly the volume fraction of solid material is 
extremely small, of order 10~ 6 . Due to this low volume 
fraction, dusty plasma is one of the softest substances 
known; judged by its shear modulus when in crystalline 
form [15], it is nineteen orders of magnitude softer than 
metals, and six orders of magnitude softer than colloidal 
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crystals [TU]. Because of its extreme softness, a crys- 
talline dusty plasma has a sound speed of only a few 
mm/s [21], and it is possible to generate velocity gradi- 
ents as large as the thermal velocity divided by the in- 
terparticle spacing, as we will show here. While a dusty 
plasma has a very small viscosity 77 [23 , it also has a 
very low mass density p, so that the kinematic viscos- 
ity v = rj/p is actually comparable to that of water, 
1 mm 2 /s [25] . 

With these extreme properties, it is worth asking 
whether significant viscous heating and a large value of 
Br can be attained in dusty plasma, and in this Letter we 
find that this is so. Our experiment is designed to yield 
all the parameters required to determine Br, Pr, and Ec. 
Our measurements include a simultaneous quantification 
of the velocity profile v x (y), a profile of a kinetic tem- 
perature, and values for two transport coefficients, r\ and 
K. 

Our experimental setup was similar to the one de- 
scribed in |22j for another experiment. A vacuum 
chamber was filled with 15.5 mTorr argon gas at room 
temperature, which was partially ionized by applying 
214 V peak-to-peak voltage at 13.56 MHz. About 10 4 
dust particles, which were 8.09 fan diameter melamine- 
formaldehyde microspheres with a mass of m<j = 4.18 x 
10~ 13 kg, were introduced into the plasma. They gained 
a large negative charge of —9700 e where e is the elemen- 
tary charge as determined by a wave-spectra method [2"5] 
which has an accuracy of about ±10%. They became 
electrically levitated by a vertical dc electric field above 
a horizontal electrode due to its dc self-bias of —138 V. A 
weaker radial dc electric field in the plasma prevented the 
charged dust particles from escaping in the horizontal di- 
rection. The dust particles filled only a single horizontal 
layer, and significant out-of-plane displacements |26j were 
not present, so that particles within the dust layer can 
be described as having mainly two-dimensional (2D) mo- 
tion. As a single dust particle moved at a velocity V^, it 
experienced a frictional drag force — ^gas^dV^ due to col- 
lisions with the gas atoms [57j, where v Eas — 2.7 s _1 . The 
dust particles had an interparticle distance b — 0.50 mm, 
corresponding to a Wigner-Seitz radius [55] a = 0.26 mm. 
At this interparticle distance, the longitudinal sound 
speed 121] in the crystalline lattice was measured to be 
ss 25 mm/s. 

We generated a stable subsonic laminar flow pattern 
analogous to a Couette flow. In a traditional Couette 
flow, which fills a 3D volume, a pair of planar plates 
manipulate a fluid between them by moving at a steady 
speed in the x direction, so that the flow velocity v has 
only a single vector component, v x (y), which varies only 
with the coordinate y perpendicular to the flow [9]. The 
shear dv x /dy is responsible for the viscous heating in 
Eq. g). 

In our experiment, instead of planar plates we used 
laser manipulation [231 123 I29"H53] of the dust particles 
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FIG. 1: (Color online) . Sketch of laser-driven flows in the 2D 
dusty plasma. In the region of interest, the flow is straight, 
with curved flow limited to the edge of the dust layer, as 
sketched in (a). A video image of the dust particles within 
the region of interest is shown in (b). In our data analysis, 
we exploit the symmetry of the experiment, and divide the 
region of interest into 89 long narrow rectangular bins, with 
one example bin shown in (b). The particle data are converted 
to continuum data by averaging data for individual particles 
in each bin. 



to cause them to flow. Laser beams can drive significant 
flows, with strong gradients, due to the extreme softness 
of the dusty plasma. Using a rapidly-rastered scanning 
mirror, we dispersed a 532-nm laser beam (with a power 
of 2.28 W as measured inside the chamber) into a sheet 
of width approximately 4 mm, which was incident on the 
dust layer at a downward angle of 6°. This laser sheet 
drove a flow in the +x direction. A second laser sheet, 
separated by about 5 mm from the first, was directed 
oppositely to drive a flow in the —x direction, as shown 
in Fig. 1. Using laser sheets of « 4 mm width (instead 
of narrowly-focused laser beams) allowed us to generate 
wider flow velocity profiles. The video can be seen in the 
Supplemental Material 34 . Our shear region allowed us 



to measure the velocity profile v x (y) and the kinetic tem- 
perature profile Tki n (y) with sufficient spatial resolution 
to characterize the three dimensionless parameters: Br, 
Pr, and Ec. 

To determine the flow velocity v and kinetic tempera- 
ture Tkin of the flowing dust particles as functions of y, 
we used video imaging and particle tracking [35 , 36J . The 
dust particles were imaged from above at 55 frames per 
second, and in each video frame we analyzed the particles 
within a region of interest measuring 23.5 mm x 23.5 mm 
and containing w 2500 dust particles. We measured [35] 
the velocity of individual dust particles. We then con- 
verted these particle data to continuum data by a spatial 
average of the particle data within rectangular bins of 
width Ay = 0.26 mm extending across the full region of 
interest [21] • In doing this, we assumed that x is an ignor- 
able coordinate in the flow, which we have validated [37] , 
and that our results for v and T kin are functions only of 
y 38 . We used the particle velocities V<; to calculate 
the flow velocity v = and the kinetic temperature 
Ikin(y) = rn<i\Vd — v| 2 /2/cs, where the overline indicates 
a spatial average. We note that the kinetic temperature, 
which is calculated using the mean-square velocity after 
subtracting the local flow velocity, is not necessarily the 
same as a thermodynamic temperature T. 

In our results for the velocity profile, Fig. 2(a), we see 
broad peaks, which are laser-driven flows. On each side of 
the broad peaks there is a region of high shear, where the 
gradient dv x /dy reaches its maximum of 4.3 s _1 . Most 
of our attention in this Letter is devoted to these shear 
regions. 

Our experiment achieved a close approximation to the 
symmetry of a Couette flow, which would ideally have 
v y /v x — everywhere and a uniform number density. 
In our experiment, \v y /v x \ f=a 10~ 2 as can be seen by 
comparing panels (a) and (b) of Fig. 2, and the number 
density in Fig. 2(c) is uniform within our measurement 
uncertainties. We note that unlike a traditional Couette 
flow for a single-phase liquid, ours has a nonlinear veloc- 
ity profile v x {y); this is due to frictional drag applied to 
the dust particles by the gas [23] . 

The first chief result of this Letter is the appearance 
of profound peaks in the kinetic temperature profile, 
Fig. 2(d). These peaks coincide with the locations of 
maximum shear, \dv x /dy\. Unlike most cases where vis- 
cous heating in a shear flow leads to only a tiny change 
in temperature, here we see a very large effect, with a 
six-fold change in absolute temperature over only a few 
mm. We will discuss the significance of these tempera- 
ture peaks below. 

In the corresponding video of the particle motion [31] , 
one can see the effect of viscous heating in the motion of 
individual particles. In the high shear region, there are 
frequent collisions of particles at slightly different values 
of y, so that they are flowing at different speeds, and 
these collisions cause a scattering of momentum and en- 
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FIG. 2: (Color online) . Spatially-resolved profiles of the con- 
tinuum quantities. The two broad peaks in the profile of v x 
in (a) are due to the laser manipulation. We are interested 
mainly in the regions of high shear, dv x /dy. The symme- 
try of the experiment minimized any flows in the y direction, 
(b), and the number density was nearly uniform, (c). One of 
our chief results is the appearance of profound temperature 
peaks, (d). This kinetic temperature, shown in Kelvin, was 
determined directly from the motion of individual dust par- 
ticles. The profiles of v x and Tki n allow us to determine Br 
and Ec, and verify that viscous heating has a strong effect. 
Marked at the bottom are the regions of high shear, where 
\dv x /dy\ > 1.3 s" 1 . 

ergy that increases the random fluctuations of particle ve- 
locities. These particle velocity fluctuations correspond 
directly to a high kinetic temperature. This is the viscous 
heating effect, as observed at the particle level. 

In order to achieve large values of Br and Ec, it is de- 
sirable to have a velocity shear that is large in dimension- 
less units. To make it dimensionless, we divide dv x /dy by 
vth/d, where Vth — ^J^bT jma. is the thermal velocity 
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and a is the the Wigner-Seitz radius. The dimensionless 
shear (dv x /dy) /(vth/a) in our experiment is 1. This is 
a large value; for comparison, the dimensionless shear is 
about 10 -4 for a laminar boundary layer of air flowing 
at 100 m/s over a cm-size flat plate [55] . 

To determine the three dimensionless numbers (Br, Pr, 
and Ec), we require not only the spatial profiles of flow 
velocity and temperature in Fig. 2, but also values for 
transport coefficients rj and k. We obtained these coeffi- 
cients by minimizing residuals of both the energy equa- 
tion, Eq. Q , and the momentum equation, while treating 
r\ and k as free parameters. The residual is an equa- 
tion's right-hand-side, which is finite (instead of zero, 
as in the ideal case) when using an input of experimen- 
tal data. Details of this residual-minimization method 
are presented elsewhere [37], where we obtained our re- 
sults for the kinematic viscosity and thermal diffusivity, 
v = r/f p = 0.69 mm 2 /s and x = k /( c p) — 8 mm 2 /s, 
respectively. These values, which were obtained simul- 
taneously from the same experiment, are consistent with 
values reported for v in a previous dusty plasma experi- 
ment [23] and \ m another [40] . 

As the second chief result of this Letter, we now de- 
termine three dimensionless numbers for the flow in our 
dusty plasma. 

The Prandtl number Pr is an intrinsic property of a 
substance. We obtain its value by evaluating Eq. ([2| 
using the values for v and x given above. We find Pr = 
v/x = (0.69 mm 2 /s)/(8 mm 2 /s) = 0.09 for our dusty 
plasma. For comparison, typical values for air |41j and 
glycerin [3] are Pr « 0.7 and 11000, respectively. 

The Eckcrt number Ec characterizes dissipation in a 
flow. Unlike Pr, it is not an intrinsic property of ma- 
terials, but depends on the flow profiles. We obtain its 
value using Eq. ^ with data from Figs. 2(a) and 2(d). 
Examining Fig. 2(a), we find Av — 2.5 mm/s, where 
we specified Av as the difference of the flow velocity be- 
tween the positions y = 8.3 mm and y = 11.0 mm, cor- 
responding to maximum shear \dv x /dy\ and minimum 
speed \v x \, respectively. Using Fig. 2(d), we obtain the 
kinetic temperature difference AT^m, specified for the 
same two positions. From Fig. 2(d) we find ATki n = 
4.6 x 10 4 K - 1.3 x 10 4 K = 3.3 x 10 4 K. We then calcu- 
late cAT ki „ = A(vf h /2) = k B AT Ma /m d = 1.1 mm 2 /s 2 . 
Inserting these values in Eq. ([3]), we obtain Ec = 5.7. 
This large value of Ec indicates that in our dusty plasma 
flow the transfer of directed kinetic energy to thermal 
energy is significant. 

The Brinkman number, which characterizes the com- 
petition between viscous heating and thermal conduction 
for a given flow, is obtained simply by multiplying the 
Prandtl and Eckert numbers, Br = PrEc. Using the val- 
ues given above, we find Br = 0.5 for the flow in our 
dusty plasma. This large value indicates that the first 
term in Eq. Q is generally half as large as the second 
term, i.e., viscous heating is half as large as thermal con- 



duction. The appearance of the temperature peaks in the 
region of sheared flow in Fig. 2(d) is consistent with the 
large value of Br. In our search of the literature mention- 
ing the Brinkman number, we found that temperature 
peaks have certainly been predicted theoretically, for ex- 
ample in Couette flows [2] , but we found no reports of an 
experiment with a spatially-resolved in-situ temperature 
measurement like ours. The microscopic mechanism for 
viscous heating is scattering of momentum by collisions 
of particles flowing at different speeds in neighboring por- 
tions of a sheared flow, which introduces random kinetic 
energy. In many substances, this local increase in ran- 
dom kinetic energy would be rapidly transported away by 
thermal conduction, so that no significant local heating 
of the substance can be observed. In our dusty plasma, 
however, temperature peaks in the region of maximum 
shear are easily observed, indicating that thermal con- 
duction is not sufficient to flatten the temperature pro- 
files. 

To demonstrate the significance of viscous heating in 
our dusty plasma, we can compare it with other sub- 
stances using the Brinkman number. For water flow- 
ing in a channel with a width of about 1 mm (which 
is comparable to the scale lengths in our experiment), 
Tso and Mahulikar 5 found Br in the range 1 x 10~ 8 to 
17 x 10~ 8 . Such low values indicate that viscous heat- 
ing causes very little temperature change, due to the 
strong competing effect of thermal conduction. For glyc- 
erin in a Taylor-Couette flow, White and Muller found 
Br = 0.0359 4J. For a polymer flowing as it is squeezed 
in an extruder used for manufacturing thermoplastics, 
where viscous heating is important because it melts the 
substance, Housz and Meijer [42] found Br = 0.47. In 
our literature search, the only experiment that we found 
with a Br larger than ours was performed using an ex- 
otic liquid that requires a different definition of Br [43] . 
We have also not found any previous experimental study 
of viscous heating that includes a measurement of the 
temperature profile, revealing the spatial localization of 
viscous heating, as we have done in Fig. 2(d). 

In summary, using in-situ measurements of the kinetic 
temperature in a highly sheared flow of a dusty plasma, 
we have demonstrated that viscous heating can result in 
profound temperature peaks. These temperature peaks 
coincide with the locations of maximum shear. We deter- 
mined Br = 0.5, Pr = 0.09, and Ec = 5.7. Our large value 
of Br indicates viscous heating to an extent that is seldom 
measured experimentally in other substances. Further 
details of the experiment, including spatially resolved de- 
terminations of the terms in the continuum equations for 
momentum and energy, will be presented in another pa- 
per [37] . 
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